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ABSTRACT 

The  functional  dependence  of  the  rate  constants  for  self  exchange,  k>x , 
for  a  series  of  metallocene  redox  couples  to  solvent -induced  variations  in  the 
nuclear  frequency  factor,  i/n,  engendered  by  alterations  in  the  longitudinal 
solvent  relaxation  time,  rL,  are  utilized  to  deduce  values  of  the  electronic 
matrix  coupling  element,  H12,  for  electron  exchange.  The  analysis  exploits 
the  sensitivity  of  the  k<x  -  r^1  dependence  to  the  degree  of  reaction 

adiabaticity  and  hence  H12  for  a  given  electron- exchange  reaction.  Six 

metallocene  couples  are  examined:  Cp2Co+/0,  Cp2Fe+/0  (Cp  -  cyclopentadienyl) , 
and  the  decamethyl  derivatives  Cp2Co+/0  and  Cp2Fe+/0  scrutinized  previously 
(refs.  5a, b) ,  additional  solvent -dependent  k#x  values  for  carboxy- 
methyl(cobaltocenium-cobaltocene)  [CpJCo+/0,  e  -  "ester"]  and  hydroxy - 
methyl(ferrocenium-ferrocene)  (HMFc+/0J.  Kinetic  data  are  examined  in  15 
solvents,  including  11  "Debye"  solvents  for  which  it  is  anticipated  that 
»/n  a  r^1 .  Corrections  to  k#x  for  the  solvent -dependent  variations  in  the 
barrier  height  were  obtained  by  corresponding  measurements  of  the  optical 
electron- transfer  energies  for  the  related  binuclear  complex 
biferrocenylacetylene .yielding  "barrier- corrected"  rate  constants,  k^x.  The 
k^x  -  r*1  dependencies,  as  well  as  the  k#x  values  in  a  given  solvent,  are 

markedly  dependent  on  the  redox  couple.  The  log  k^x  -  log  r*1  plots  for  the 
most  facile  couple,  Cp2Co+/o,  exhibit  slopes  approaching  unity  for  smaller 
f^1  values.  The  less  facile  couples  yield  smaller  slopes,  diminishing  in  the 
same  sequence  that  k^x  decreases  in  a  given  solvent:  Cp2Co+/0  >  Cp*Co+/0  > 

Cp2Co+/°  >  Cp2Fe+/0  >  Cp2Fe+/0  >  HMFc+/0.  These  findings  are  consistent  with 
H12  decreasing  in  the  same  order.  Comparison  of  such  rate -solvent  friction 
dependencies  with  corresponding  plots  calculated  using  a  combined  solvent 

friction-electron  tunneling  model  yield  the  following  approximate  values  of 
the  matrix  coupling  element  for  reactant  "closest  approach",  H°2  (kcal  mol*1): 
Cp^Co+/°,  1.0;  Cp;Co+/°,  0.5 -1.0;  Cp2Co+'°,  0.5;  Cp^Fe+'°,  0.2;  Cp2Fe+'°,  0.1; 

HMFc+/0,  0.075.  Reasonable  concordance  is  seen  with  recent  theoretical 
estimates  of  H®2  for  Cp2Co+/°  and  Cp2Fe+/0.  The  relationship  between  H°2  and 
metallocene  electronic  structure  is  briefly  discussed.  The  analysis  also 
enables  effective  solvent  relaxation  times  for  adiabatic  barrier  crossing  in 
non-Debye  media,  including  primary  alcohols,  to  be  extracted. 
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A  central  fundamental  question  in  electron  transfer  of  widespread  current 

interest  concerns  how  the  extent  of  overlap  between  the  donor  and  acceptor 

redox  orbitals,  and  its  role  in  influencing  reaction  rates,  depends  upon  the 

electronic  properties  and  spatial  configuration.  The  overlap  is  described  by 

the  electronic  matrix  coupling  element,  H12,  which  determines  the  extent  of 

1  2 

resonance  splitting  of  the  barrier  top.  ’  This  factor  exerts  a  crucial 
influence  upon  electron- transfer  kinetics,  not  only  by  influencing  the  barrier 
shape  and  height  but  also  the  degree  of  reaction  adiabaticity  (i.e.,  the 
probability  that  electron  transfer  will  occur  once  the  transition  state  has 
been  formed).  Recent  ab  initio  calculations  of  H12  for  simple  outer-sphere 
electron  exchange  processes  illustrate  the  sensitivity  of  the  coupling 
strength  to  electronic  structure  as  well  as  to  the  spatial  configuration  of 

O 

the  reacting  pair.  While  experimental  estimates  of  H12  for  optical  electron 

la 

transfer  within  binuclear  complexes  can  be  extracted  from  band  intensities, 
the  evaluation  of  the  matrix  coupling  element  for  thermal  outer- sphere 
electron  transfer  is  considerably  less  straightforward.  In  spite  of  the 

importance  of  the  latter  quantity  as  alluded  to  frequently  in  the 
literature,  its  experimental  evaluation  has  remained  largely  elusive.^ 

We  have  recently  been  examining  in  detail  the  solvent-dependent  electron- 
exchange  kinetics  of  various  metallocene  redox  couples  in  homogeneous 
solution^  as  well  as  electrode  surfaces^  in  order  to  explore  the  role  of 
solvent  relaxation  dynamics  on  the  barrier-crossing  frequencies.  Of  primary 
concern  have  been  redox  couples  of  the  form  Cp2M+/0  or  Cp2M+/o,  where  M  -  Co 
or  Fe,  Cp  -  cyclopentadienyl ,  and  Cp'  -  pentamethylcyclopentadienyl , ^ ^ 
although  other  metallocene  couples  have  also  been  examined. During  these 
studies,  it  became  apparent  that  the  rate  constants  for  homogeneous  self 
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exchange,  k#Jt,  in  a  given  solvent  display  a  marked  sensitivity  to  the 

metallocene  electronic  structure.  In  particular,  k#Jt  values  for 

cobaltocenium-cobaltocene  couples  are  substantially  (up  to  ca  10  fold)  larger 

than  for  corresponding  ferrocenium- ferrocene  couples:  up  to  100  fold 

enhancements  in  k<Jt  are  observed  upon  methylation  of  the  Cp  rings  together 

5a  7 

with  substituting  cobalt  for  iron.  *  These  rate  variations,  especially 
those  resulting  from  metal  substitution,  were  traced  primarily  to  the  greater 
extent  of  donor-acceptor  orbital  overlap  for  the  cobalt  versus  iron  systems 
anticipated  from  the  primarily  ligand-  versus  metal-centered  character  of  the 
redox  orbitals.^  These  differences  in  orbital  overlap  and  hence  H12  are 
consistent  both  with  optical  electron- transfer  measurements  for  bicobaltocene 

Q 

versus  biferrocene  cations  and  with  recent  ab  initio  calculations  (vide 
infra) . ® 

Generally  speaking,  the  magnitude  of  H12  for  a  particular  internuclear 
configuration  (i.e.,  precursor  complex  geometry)  can  affect  the  rate  constant 
for  that  geometry,  k<t  (s"1),  by  influencing  the  electronic  transmission 
coefficient,  /c#1,  and  additionally  the  free-energy  barrier,  AG*,  as  in  the 
conventional  expression^ 

k.t  "  *.i  e*p(-AG*/RT)  (1) 

where  vn  is  the  nuclear  frequency  factor.  The  last  term  describes  the  net 
dynamics  along  the  nuclear  reaction  coordinate  in  the  vicinity  of  the  barrier 
top,  which  may  include  contributions  from  both  collective  solvent  and  inner - 
shell  motion.  If  only  a  narrow  range  of  precursor  geometries  contributes  to 
the  bimolecular  reaction  (vide  infra),  k#t  can  be  related  simply  to  kfx  by 


(2) 
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where  Kj,  is  the  statistical  probability  of  forming  the  precursor 
geometry .^’10  If  HJ2  is  very  small,  the  reaction  will  be  nonadiabatic  such 

It  O 

that  «  1  and  k-x  a  H22.  ’  In  the  adiabatic  regime  where  H12  is 

sufficiently  large  so  that  -*  1,  the  dependence  of  k>t  (and  hence  k<x)  upon 

H12  will  be  milder,  arising  primarily  from  the  diminution  in  the  barrier 

Xc  3  4&  o 

height.  Evidence  from  both  experimental  ’  ’  and  theoretical*  sources 
suggests  that  bimolecular  reactions ,  even  involving  small  inorganic  reactants , 
are  often  nonadiabatic  in  nature . 

While  estimates  of  and  hence  H12  could  in  principle  be  obtained  from 

k#x  measurements  for  a  single  system  by  using  Eqs.  (1)  and  (2),  uncertainties 

in  the  required  estimates  of  AG*  and  Kp  make  this  procedure  of  little  utility. 

One  approach  is  to  examine  the  relative  k#x  values  for  a  series  of 

structurally  related  reactions  chosen  so  that  differences  in  the  extent  of 

electronic  coupling  provide  a  likely  major  contribution  to  the  observed  rate 
4a  7 

variations.  ’  However,  due  to  the  involved  dependence  of  /c#1  upon  H12  this 

procedure  requires  at  least  an  estimate  of  k#x  when  x#1  -»  1  (vide  infra) . 

An  attractive,  yet  virtually  unexplored,  possibility  is  to  examine 

additionally  the  sensitivity  of  k#x  to  variations  in  the  nuclear  frequency 

factor  This  is  because  the  dependence  of  k>x  upon  will  be  sensitive  to 

the  degree  of  reaction  nonadiabaticity ,  which  in  turn  is  influenced  by  the 

11  12 

magnitude  of  as  well  as  H12.  ’  Thus  while  we  expect  that  kBX  a  i/n  under 

adiabatic  conditions  [i.e.,  when  -  1,  Eq.  (1)],  increasing  i/D  for  a  given 
H12  eventually  can  lead  to  a  nonadiabatic  condition  whereupon  k<x  becomes 
Independent  of  vn,  since  then  «cxl  a  v*1.^  The  form  of  the  ktx-i/n  dependence 
should  therefore  be  sensitive  to,  and  therefore  diagnostic  of,  the  magnitude 
of  H12.^  While  i/n  for  processes  featuring  large  inner- shell  distortions  will 
be  determined  at  least  partly  by  reactant  vibrations  that  cannot  be  varied 
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readily,  the  dynamics  of  collective  solvent  motion  can  provide  the  dominant 

contribution  to  when  the  inner- shell  barrier  is  small  (say  <  1  kcal 
13  *15 

mol'1).  *  Substantial  (up  to  20-50  fold)  variations  in  i/n  can  readily  be 

induced  under  the  latter  conditions  by  judicious  alterations  in  the  solvent 
medium,  since  in  many  cases  i/n  a  r^1,  where  rL  is  tha  longitudinal  solvent 
relaxation  time.^a’^a  c’^*^  This  enables  the  dependence  of  k#Jt  upon  i/n  to 
be  examined  provided  that  the  rates  are  corrected  for  the  corresponding 
variations  in  AG*"5’^  (vide  infra). 

In  the  present  paper,  we  exploit  such  a  solvent-dependent  analysis  to 
obtain  H12  estimates  for  a  series  of  six  ferrocene  and  cobaltocene  self- 
exchange  reactions.  The  kex  values  were  obtained  using  the  proton  NMR  line- 
broadening  technique,  with  the  required  solvent  dependence  of  AG  being 
extracted  from  optical  electron- transfer  energies  for  a  related  binuclear 
system,  biferrocenylacetylene  cation  (cf  ref.  5).  Some  of  the  data  used  in 

5 

the  present  analysis  has  been  discussed  previously  in  a  related  context. 
However,  the  acquisition  of  additional  kinetic  and  optical  barrier  data  along 
with  our  recent  reevaluation  of  the  solvent-dependent  kinetics  for 

ferrocenium- ferrocene  self  exchange'**5  enable  the  role  of  orbital  overlap  in 
electron  transfer  to  be  addressed  here  quantitatively  for  the  first  time. 
Besides  the  parent  metallocene  couples  Cp2Fe+/0  and  Cp2Co+/°,  and  the 
decamethyl  derivatives  Cp2Fe+/0  and  Cp2Co+/°,  we  utilize  new  solvent -dependent 
rate  data  for  the  carboxymethyl(cobaltocenium-cobaltocene)  and  hydroxy- 
methyl(ferrocenium-ferrocene)  couples,  (CpC02Me)2Co+/o  and  (Cp ■  CpCH20H)Fe+/o , 
respectively , '5C  abbreviated  here  as  "Cp2Co+/°"  (e  -  ester)  and  "HMFc+/c". 

These  latter  two  systems  were  selected  initially  in  view  of  their  water 
solubility,  enabling  barrier-crossing  dynamics  in  aqueous  and  nonaqueous  media 
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to  be  compared.  Together,  these  six  reactant  systems  are  shown  to  span  a 
wide  (ca  tenfold)  range  in  H12,  from  ca  0.1  to  1.0  kcal  mol"1. 

EXPERIMENTAL  SECTION 

Carboxymethylcobaltocenium  tetrafluoroborate  was  synthesized  by  stirring 
the  corresponding  carboxylic  acid  derivative  (CpC00H)2Co‘BF4) ,  prepared  as  in 
ref.  16,  in  methanol  purged  with  anhydrous  HC1  for  24  hours.  Addition  of  HBF4 
and  evaporating  some  solvent  precipitated  the  desired  product,  which  was 
recrystallized  from  hot  methanol.  The  reduced  species, 

carboxymethylcobaltocene ,  was  prepared  by  stirring  a  hexane  suspension  of  the 
salt  with  a  deficiency  of  cobaltocene.  The  solution  was  filtered  to  remove 
Cp2Co*BF4  and  the  filtrate  evaporated  to  dryness  under  vacuum. 
Hydroxymethylferrocene  was  synthesized  from  trimethylaminoferrocene  (Aldrich) 
as  described  in  ref.  17.  The  oxidized  form,  (Cp • CpCH2OH) Fe ■ PF6 ,  was  prepared 
by  shaking  a  hexane  solution  of  the  reduced  species  with  an  acidic  aqueous 
solution  containing  a  twofold  excess  of  FeCl3  until  the  hexane  phase  became 
colorless.  After  filtering,  a  solution  of  NH4PF6  was  added  to  the  aqueous 
phase  to  precipitate  the  desired  product.  Biferrocenylacetylene  (BFA)  was 
prepared  as  described  in  ref.  18;  the  acetylferrocene  precursor  was  obtained 
by  acetylation  of  ferrocene  as  outlined  in  ref.  19.  For  all  complexes, 
identification  and  purity  were  confirmed  by  proton  NMR  and  cyclic  voltammetry. 

The  near-infrared  spectra,  so  to  obtain  optical  electron- transfer 

on  oo 

energies  for  biferrocenylacetylene  cation  (BFA+) ,  ’  were  measured  between 

800  and  2000  nm  with  a  Cary  Model  17D  spectrophotometer,  using  1  cm  quartz 
cells.  Typically,  ca  1  mil  concentrations  were  employed  so  to  yield  maximum 
absorbances  of  about  0.2-0. 5.  For  most  solvents,  BFA+  was  generated  from  BFA 
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by  in-situ  oxidation  with  one  equivalent  of  Fe(bpy)3(PF6)3  (bpy  - 
2 , 2 'bypyridyl) .  In  some  cases  (e.g.  methanol),  the  solid  BFA+BF4  salt  was 
employed.  This  was  isolated  by  oxidizing  BFA  in  nitromethane  under  nitrogen 
with  a  stoichiometric  quantity  of  AgBF4  (Aldrich).  After  filtering,  the 
product  was  isolated  by  adding  diethyl  ether. 

Most  details  of  the  NMR  sample  preparation  measurements  and  line- 
broadening  data  analysis  are  given  in  refs.  5a  and  b.  Solvents  were  high- 
purity  grades  from  Burdick  and  Jackson,  Aldrich,  or  Fluka  and  were  purified 
where  appropriate  by  using  standard  procedures.  Deuterated  solvents  were 
employed  whenever  feasible.  However,  protested  media  were  also  used 
successfully  by  suppressing  the  solvent  peak(s)  using  homonuclear  irradiation 
at  the  desired  resonance  frequencies,  with  a  minimum  of  radio  -  frequency  power 
to  avoid  significant  dielectric  heating.  Proton  NMR  spectra  for  the  HMFc+/0 
and  Cp*Co+/°  systems  were  collected  on  Nicolet  NT  200  and  NT  470  instruments, 
respectively  (operated  at  200.0  and  469.5  MHz).  Due  to  peak  splitting  on  the 
Cp  ring  protons,  the  line -broadening  measurements  for  these  two  couples 
utilized  the  -CH2-  and  methyl  substituent  proton  resonances,  respectively.  A 
complication  for  the  former  is  that  the  methylene  resonance  appears  as  a 
doublet  in  the  pure  diamagnetic  species  (i.e.  hexamethylferrocene)  due  to 
spin-spin  coupling  with  the  hydroxyl  proton.  This  difficulty,  however,  was 
circumvented  by  collapsing  the  doublet  to  a  single  peak  by  saturating  the  -OH 
resonance  during  NMR  data  acquisition.  In  any  case,  essentially  equivalent 
line  widths  were  obtained  irrespective  of  -OH  peak  saturation,  presumably  due 
to  weak  CH2-0H  coupling.  A  difficulty  faced  with  Cp2Co+/0  is  that  the  contact 
shift  for  the  methyl  protons  is  relatively  small  even  at  470  MHz  field,  ca 
12  ppm,  which  combined  with  the  relatively  large  k  values  yields  relatively 
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small  electron- exchange  line  broadening.  Consequently,  it  was  advisable  to 
correct  the  measured  lin*'  widths  for  magnetic  field  inhomogeneity.  For 
nonaqueous  solvents,  this  was  achieved  simply  by  subtracting  the  corresponding 
measured  linewidth  for  the  tetramethylsilane  reference  signal;  the  correction 
was  estimated  to  be  1-2  Hz  in  D20.  This  complication  increases  the 
uncertainties  in  the  derived  k>x  values  for  the  CpJCo+/0  couple,  but  they  are 
still  probably  reliable  to  within  +  20-30%.  For  the  other  systems,  k#x  is 
typically  reproducible  to  10-20%.  Experimental  NMR  parameters  for  the 
additional  reactions  studied  here  are  summarized  in  the  Supplementary 
Material . 

RESULTS  AND  DISCUSSION 

A  summary  of  the  k<x  values  used  in  the  present  analysis,  for  six  redox 
couples  in  fifteen  solvents,  is  given  in  Table  I.  Besides  the  additional  data 
gathered  in  the  present  work,  this  summary  contains  pertinent  data  culled  from 
refs.  5a-c  (see  footnote  to  Table  I  for  specific  details).  For  convenience, 
the  k#x  values  in  Table  I  all  refer  to  a  uniformly  low  ionic  strength 
in  ~  0.01-0.02  H)  .  This  entailed  adjusting  the  values  for  Cp20o+/o  and 
Cp2Co+/°  given  earlier  in  ref.  5a,  since  these  were  obtained  at  /i  *  0.14. 
Fortunately,  this  correction  is  small  (ca  25%  in  k#x)  since  kxx  for  these 
metallocene  couples  are  insensitive  to  variations  in  ionic  strength,  at  least 
when  "noncoordinating"  anions  (such  as  PFj,  BF^)  are  employed.^  Kinetic  data 
in  dichlorome thane  and  dichloroethane ,  reported  in  ref.  5a,  are  not  included 
here.  This  is  partly  in  view  of  our  subsequent  observation  of  deviations  from 
second-order  kinetics  in  these  solvents,  probably  associated  with 
complications  from  ion  pairing.  Large  ion-pairing  effects  upon  the  reaction 
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energetics  in  such  weakly  polar  media  are  also  evident  in  optical  electron- 

transfer  data.^a  The  k#x  values  for  Cp2Fe+/0,  taken  in  part  from  ref.  5b, 

36 

differ  substantially  (up  to  ca  10  fold)  from  those  reported  by  Wahl  et  al 

which  were  utilized  in  our  earlier  solvent  dynamical  analyses. The 

origins  of  the  systematic  solvent-dependent  errors  in  these  older  data  are 

discussed  in  ref.  5b.  Experimental  difficulties,  such  as  reactant  solubility 

and  decomposition,  along  with  solvent  interferences  and  signal-to-noise 

constraints  in  the  NMR  line -broadening  measurements,  precluded  the  evaluation 

of  k#x  for  a  number  of  redox  couple/solvent  combinations.  For  example,  the 

ferrocene  couples  could  not  be  examined  in  amide  solvents  due  to 

decomposition.  Nevertheless,  sufficient  kinetic,  along  with  optical  electron 

transfer  (vide  infra) ,  data  have  now  been  assembled  to  enable  a  relatively 

detailed  dynamical  analysis  to  be  undertaken. 

Of  the  fifteen  solvents  in  Table  I,  the  first  eleven  are  considered  to 

approximate  "Debye"  dynamical  behavior  in  that  the  dielectric  loss  spectra  can 

37 

be  described  in  terms  of  a  single  relaxation  time,  rD.  In  this  case,  the 
adiabatic  barrier-crossing  frequency  i/n  for  electron-transfer  reactions  having 

little  or  no  inner-shell  barrier  is  anticipated  to  be  determined  by  the 

13  14  38  39 

longitudinal  relaxation  time,  ’  given  approximately  by  ’ 

Ti  "  <««>A0>  td  O) 

where  t.  and  «o  are  the  "infinite-"  and  zero- frequency  solvent  dielectric 
constants,  respectively.  The  values  for  each  of  these  solvents  are  given 
in  Table  I  (see  footnotes  and  refs.  5a  and  6c  for  data  sources);  the  Debye 
solvents  are  listed  in  order  of  increasing  rL.  Rate  data  obtained  in  three 
likely  "non-Debye"  solvents  -  propylene  carbonate,  methanol,  and  ethanol  -  are 
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also  Included  in  this  compilation,  along  with  r^1  estimates  for  the  major 

(low-frequency)  relaxation,  given  in  parentheses.  (Such  non-Debye  behavior  is 

characterized  by  at  least  one  higher -frequency  dispersion  region  having 

25  33-35  37 

dissipative  characteristics.  *  ’  )  Although  the  connection  between  the 

dielectric  loss  spectra  and  vB  is  considerably  less  straightforward  for  non- 
Debye  solvents, <25,40  tjjese  kiltie  data  are  nevertheless  of  interest  to  the 
present  discussion  (vide  infra) . 

As  already  noted,  the  central  thrust  of  the  present  work  is  to  examine 
the  sensitivity  of  the  barrier-crossing  frequency  to  alterations  in  the 
solvent  dynamics  as  a  means  of  probing  the  orbital -overlap  factor,  H12.  Even 
though  only  relative,  rather  than  absolute,  preexponential  factors  for  a  given 
reaction  in  a  suitable  series  of  solvents  are  required  for  this  analysis,  the 
solvent -dependent  k#x  values  still  need  to  be  corrected  for  differences  in  the 
barrier  height,  AG*.  The  required  solvent -dependent  AG*  estimates  can  be 
obtained  in  two  ways.  The  first  involves  employing  theoretical  models,  most 

41 

simply  the  dielectric  continuum  treatment  of  Marcus.  We  have  used  this 

approach  in  several  previous  discussions . ^a’ C 

An  alternative,  and  probably  more  reliable,  procedure  utilizes 

experimental  energies  for  optical  electron  transfer,  Eop,  within  symmetrical 

mixed-valence  compounds  that  are  structurally  similar  to  the  thermal  electron- 

exchange  reactions  of  interest.^  While  suitable  mixed-valence  systems  are  not 

abundant,  the  biferrocenylacetylene  cation  (BFA+)  constitutes  an  almost  ideal 

choice  for  the  present  purposes.  Besides  being  structurally  closely  similar 

to  the  reactants  of  interest  here,  the  juxtaposition  of  the  Cp2Fe+/Cp2Fe 

20 

partners  in  the  binuclear  complex  roughly  mimics  that  anticipated  for  the 
precursor  complex  for  the  thermal  self-exchange  reactions.  Most  importantly, 
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BFA+  conforms  approximately  to  "class  II"  (i.e.  valence- trapped) 

20  42  * 

behavior,  ’  enabling  solvent -dependent  thermal  barriers,  AG£ ,  to  be 

obtained  simply  from^6 

*G*  "  Eop/4  (4) 

It  should  be  noted  that  AG*  obtained  in  this  manner  is  a  "cusp"  barrier,  i.e. 

in  the  absence  of  barrier-top  roundedness  induced  by  donor -acceptor  orbital 

la  43  * 

overlap  (H12  -*  0) .  ’  The  actual  barriers ,  AG  ,  for  the  present  thermal 

reactions  will  depend  upon  the  internuclear  configuration  as  well  as  the 
metallocene  electronic  structure  (vide  infra).  However,  at  least  the  solvent 
dependence  of  AG  for  a  given  reaction,  as  required  for  the  present  analysis, 

jL 

should  follow  approximately  the  corresponding  variations  in  AG  . 

The  AG*  estimates  obtained  from  the  experimental  Eop  values  for  BFA+  in 
each  solvent  are  also  listed  in  Table  I .  Values  of  E  for  six  of  these 

op 

20 

solvents  were  reported  earlier;  the  present  values  agree  uniformly  within 

the  experimental  uncertainty  (2-3%).  Although  substantial  ion-pairing  effects 

21 

upon  Eop  have  been  observed  for  BFA+,  these  effects  are  small  or  negligible 

at  the  low  anion  concentrations  (<  5  mM)  and  relatively  polar  media  employed 

2  lb 

for  the  present  Eop  measurements.  (A  brief  comparison  between  the  ionic 

strength  dependence  of  k#Jt  and  Eop  is  provided  in  ref.  5b.)  The  AG*  value 
given  for  D20  was  estimated  from  Eop  data  for  the  biferrocene  cation,  as 
described  in  ref.  5c,  in  view  of  the  insolubility  of  BFA+  in  this  solvent. 
Values  of  Eop  could  not  be  obtained  for  five  solvents  of  interest  here 
(dimethylsulfoxide  and  the  amides)  due  to  BFA+  decomposition.  Estimates  of 
AGe  are  nevertheless  given  in  parenthesis  in  Table  I.  These  were  obtained  by 

taking  advantage  of  the  approximate  correlation  between  Eop  and  (e‘p  -  €~x), 

20  22 

where  «op  is  the  solvent  optical  dielectric  constant.  ’ 


Although  there  are 
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significant  deviations  from  the  direct  proportionality  between  Eop  and 

(t*1  -  c*1)  predicted  from  the  dielectric  continuum  model,  the  Eop  values 

estimated  in  this  manner  are  likely  reliable  to  within  at  least  ±  0.2  kcal 

mol*1.  (The  significance  of  such  comparisons  to  the  issue  of  noncontinuum 

*  22 

solvent  effects  upon  AGC  will  be  addressed  elsewhere.  )  For  simplicity, 

these  AGC  values  are  utilized  for  each  redox  couple  considered  here.  Slightly 

(up  to  25%  smaller)  values  are  anticipated  for  the  decamethyl  derivatives  in 

41 

view  of  their  larger  size  on  the  basis  of  the  usual  continuum  treatment. 
However,  given  that  only  the  solvent  dependence  of  AG*  is  important  for  the 
present  purposes,  and  that  the  span  of  values  in  Table  I  is  relatively  small 
(<  0.9  kcal  mol'1),  such  corrections  exert  only  a  small  or  negligible 
influence  upon  the  data  analysis.  Zn  any  case,  they  are  of  marginal  validity, 
especially  given  the  uncertainties  in  the  precise  precursor  geometries  that 
contribute  to  k>x  (vide  infra) . 

Rate- Solvent  Friction  Dependencies 

In  our  previous  analyses  of  solvent  dynamical  effects  for  metallocene 
self-exchange  reactions, ^  we  obtained  solvent -dependent  barrier-crossing 
frequencies,  <c#1vn,  by  inserting  the  experimental  rate  constants  into  Eqs.  (1) 
and  (2)  together  with  the  corresponding  AG*  values  and  a  suitable  estimate  of 
Kp.  The  manner  and  extent  of  the  solvent  dynamical  influence  upon  k#x  can 
then  be  ascertained  by  examining  the  dependence  of  these  inferred  values 

upon  known  dynamical  parameters,  most  simply  r^1.^ 

A  somewhat  different  analysis  is  employed  here.  While  this  simple 
procedure  is  desirable  for  some  purposes,  as  noted  above  the  present  objective 
hinges  on  examining  the  functional  sensitivity  of  the  net  barrier-crossing 
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dynamics  to  the  solvent  friction,  as  reflected  in  r^1.  A  complication  in 

employing  the  preequilibrium  treatment  [Eq.  (2)]  for  this  purpose  is  that  the 

effective  value  for  a  given  reaction  will  tend  to  vary  somewhat  with  the 

friction  rather  than  remain  constant  as  is  assumed  in  this  analysis. This 

is  because  the  distribution  of  precursor* complex  geometries,  most  simply  the 

range  of  donor -acceptor  spatial  separations  (r) ,  contributing  significantly  to 

the  bimolecular  reaction  rate  will  inexorably  alter  as  the  barrier-crossing 

dynamics  and  hence  degree  of  reaction  adiabaticity  changes. ^^,12  strictly 

speaking,  the  bimolecular  rate  constant,  k#x,  even  for  a  pair  of  "ideal" 

spherical  reactants  is  composed  of  an  integral  of  "local"  unimolecular  rate 

44 

constants,  k#t(r),  expressed  as 

on 

k„  -  <4*N/103)  /  r2  k,t(r)  gDA(r)  dr  (5) 

ro 

where  N  is  Avogadro's  number,  ro  is  the  donor -acceptor  distance  corresponding 
to  the  reactant's  closest  approach,  and  gDA(r)  is  the  radial  pair  distribution 
function.  [For  simplicity,  gDA(r)  will  be  taken  here  as  unity.]  The  radial 
dependence  of  k#t,  and  hence  k#x,  is  determined  primarily  by  two  factors:  the 
diminution  in  H12  (and  hence  of  «#1)  and  the  increases  in  AG*  expected  with 
increasing  r.  On  the  basis  of  Eq.  (1),  both  of  these  effects  lead  to  k#t(r) 
decreasing  sharply  with  increasing  r.  In  addition,  in  the  presence  of  solvent 
friction  is  anticipated  to  increase  significantly  both  with  decreasing  H12 
and  increasing  AG*,  and  hence  with  increasing  r.^’^  All  three  terms  that 
comprise  k#t  in  Eq.  (1)  are  therefore  dependent  on  the  precursor-complex 
geometry  and  can  therefore  Influence  the  k<x  -  r’1  dependence  in  a  relatively 
complex  fashion. 

Although  this  situation  may  seem  somewhat  daunting,  these  complications 
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nay  be  circumvented  in  the  following  manner.  Since  the  AGC  values  in  Table  I 
span  the  narrow  range  4.55  to  5.4  kcal  mol*1,  the  solvent  dependence  of  k#x 
due  to  variations  in  the  barrier  height  may  be  corrected  for  by 

Kx  -  Kx  **pKAG*  -  5 . 0)/RT]  (6) 

where  R  is  the  gas  constant  and  k^s  is  the  rate  constant  that  would  be 

observed  in  a  given  solvent  if  the  cusp-limit  barrier  was  5.0  kcal  mol*1. 

Since  the  deviations  from  this  average  AGe  value  for  the  present  solvents  are 

less  than  0.5  kcal  mol*1,  the  corrections  to  k#x  amount  to  only  twofold  or 

less.  Although  the  inevitable  uncertainties  in  AG*,  perhaps  as  much  as  20-25% 

(vide  supra) ,  in  the  actual  AGC  values  may  engender  some  error  in  the  derived 
45 

k^x  values,  at  least  the  required  relative  kxx  values  in  the  various 
solvents  should  be  reliable  to  within  50%  or  less.  The  solvent  dependence  of 
kxx  will  therefore  reflect  chiefly  variations  in  the  preexponential  factor 
caused  by  alterations  in  the  solvent  friction.  The  nature  and  extent  of  this 
dependence  in  turn  should  be  sensitive  to,  and  hence  diagnostic  of,  the  degree 
of  reaction  adiabaticity ,  and  hence  Hu,  for  the  precursor-complex 
configurations  that  contribute  most  prominently  to  the  observed  reaction 
rates.  This  procedure  is  related  closely  to  our  previous  analysis,^3 
involving  the  estimation  of  solvent -dependent  Kulvn  values  from  the  measured 
rate  constants  by  employing  the  preequilibrium  treatment  [Eqs.  ( 1 ) , < 2 ) ]  and 
correcting  for  the  barrier  height.  Such  estimates  are  instructive  for 

some  purposes. ^  However,  the  spatial  integration  embodied  in  Eq.  (5) 
necessarily  considers  a  range  of  precursor  geometries  with  spatially  dependent 
preexponential  and  exponential  factors  (vide  infra) .  For  the  present 
purposes,  therefore,  it  is  preferable  to  explore  the  solvent-dependent 
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dynamical  effects  in  terms  of  k^x  rather  than  to  extract  t/n  estimates  from  the 
rate  data. 

Plots  of  log  k'x  versus  log  for  five  metallocene  couples  in  eleven 
■Debye"  solvents,  extracted  from  the  k#x  and  AG*  values  in  Table  I  by  using 
Eq.  (6),  are  given  in  Fig.  1.  The  solvent  numbering  scheme  is  as  listed  in 
Table  I;  the  three  filled  symbols  (circles,  squares,  and  triangles)  refer  to 
Cp^Co+/°,  Cp2Co+/°,  and  Cp2Co+/°,  respectively,  whereas  the  two  open  symbols 
(triangles  and  squares)  refer  to  Cp2Fe+/o  and  HMFc+/°,  respectively.  (The 
points  for  Cp2Fe+/o  are  not  shown  for  clarity  in  Fig.  1  since  these  are 
uniformly  close  to  the  corresponding  Cp2Co+/o  points.)  Comparably  shaped 
log  k^x  -  log  r*1  plots,  yet  having  somewhat  larger  shapes,  are  also  obtained 
if  the  dielectric  continuum  estimates  of  AGC  (as  used  in  ref.  5a)  are  employed 
instead  to  estimate  k^x  from  k#x.  As  noted  above,  however,  the  present 
experimentally  derived  values  of  AGC  are  probably  more  reliable  than  such 
theoretical  estimates. 

Inspection  of  Fig.  1  reveals  that  the  degree  of  dependence  of  log  k^x 
upon  log  r'1  depends  markedly  upon  the  redox  couple.  For  the  most  facile 
Cp2Co+/°  couple,  the  log  k^x  -  log  r^1  slope  approaches  unity,  yet  appears  to 
decrease  as  the  solvent  friction  decreases  (i.e.  for  higher  r^1  values).  A 
similar  log  k^x  -  log  r^1  dependence  is  obtained  for  CpJCo+/t>,  although  the 
data  set  is  more  restricted.  The  less  facile  Cp2Co+/o  couple  also  displays  a 
qualitatively  similar  dependence  of  log  kxx  upon  log  r^1;  however,  the  slopes 
are  noticeably  smaller  (Fig.  1).  While  the  k^x  value  for  Cp2Co+/0  in 
benzonitrile  (point  8)  is  significantly  (ca  twofold)  larger  than  expected  on 
this  basis,  the  remaining  log  k^x  -  log  r'1  correlation  is  satisfactory, 
displaying  a  decidedly  fractional  slope  (ca  0.3-0. 5)  towards  the  highest  r*1 
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values.  The  least  facile  ferrocene  couples  Cp2Fe+/°  and  HMFc+/0  exhibit  no 
clear  dependence  of  log  k^x  upon  log  r'1 ,  even  though  the  absence  of  rate  data 
for  these  reactions  in  amides  restricts  somewhat  the  range  of  solvents. 

Consequently,  therefore,  an  interesting  spectrum  of  dynamical  behavior  is 
evident.  This  spans  the  range  from  apparently  full  "solvent  dynamical 
control",  where  kxx  a  r"1,  for  the  most  facile  cobaltocene  reactions  in 
relatively  high  friction  (small  r^1)  media,  to  largely  nonadiabatic  behavior 
("electron- tunneling  control"),  reflected  in  the  insensitivity  of  k^x  to  r'1 
for  the  least  facile  ferrocene  systems.  As  outlined  below,  the  combination  of 
such  varying  rate-friction  dependencies  for  such  a  "homologous  series"  of 
reactions,  together  with  the  rate  differences  between  the  reactions  in  a  given 
solvent,  enable  absolute  estimates  of  Hl2  to  be  obtained.  Uncertainties  in 
the  data  analysis  together  with  other  factors  limit  our  confidence  in  some 
quantitative  details  (vide  infra).  Nevertheless,  the  solid  traces  included  in 
Fig.  1  represent  intuitive  "best  fits"  of  the  log  k^x  -  log  r*1  dependencies 
for  each  reaction.  (Only  one  trace  is  drawn  for  the  Cp2Fe+/0  and  HMFc+/° 
couples  in  view  of  their  closely  similar  kinetic  behavior). 

These  results  are  displayed  in  a  somewhat  different  form  in  Table  II, 
which  lists  k^x  values  ratioed  to  the  corresponding  quantity  in  acetonitrile, 
k^x(ACN).  (Acetonitrile  was  chosen  as  a  "reference"  solvent  here  in  view  of 
its  rapid  relaxation  behavior  together  with  the  availability  of  reliable  k#x 
values  for  all  six  metallocene  couples.)  Given  alongside  these  k^x/k^x(ACN) 
values  are  the  corresponding  r^/r^1  (ACS')  ratios,  which  should  approximate  the 
rate  ratios  only  if  overdamped  solvent  dynamical  control  is  maintained.  While 
the  k^A^CACN)  values  for  Cp^Co*'0  and  Cp2Co+/o  self  exchange  in  Debye 
solvents  decrease  uniformly  in  parallel  with  r^/r^1  (ACN) ,  the  rate  ratios  for 
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the  other  couples  in  a  given  solvent  are  seen  to  increase  in  the  sequence 

Cp^Co+/0  -  CpJCo*'0  <  Cp2Co+/0  <  Cp^Fe+/°  <  Cp2Fe+/°  <  HMFc+/0,  the  dependence 

becoming  more  marked  as  t£1/t^1(ACN)  decreases  (Table  II).  This  finding  is 

consistent  with  a  progressive  decrease  in  the  degree  of  solvent  dynamical 

control  for  the  less  facile  redox  couples,  the  dependence  of  k^x  upon  r'1 

being  almost  entirely  absent  for  Cp2Fe+/°  and  HMFc+/0  (cf  ref.  5b). 

Also  included  in  Table  II  are  corresponding  data  for  the  three  non- Debye 

solvents  propylene  carbonate  (PC),  methanol,  and  ethanol,  along  with  the 

strongly  hydrogen-bonded  liquid  N-methylformamide  (NMF) .  In  accordance  with 

5a  25 

previous  discussions,  ’  the  k^x  values  for  Cp2Co+/0  and  Cp2Co+/0  in  these 

solvents  are  much  largei  than  expected  from  the  corresponding  r'1  values.  By 

comparing  the  k^x  values  with  corresponding  values  in  Debye  media  so  to 

achieve  the  best  consistency  with  the  other  k^-r^1  data,  "effective"  inverse 

relaxation  times ,  r.ff  can  be  deduced  in  these  non-Debye  solvents  (cf 

ref.  40d) .  This  procedure  yields  approximate  values  as  follows:  PC, 

2  x  1012  s'1;  methanol,  >  3  x  1012  s'1;  ethanol,  4  x  1012  s'1;  NMF,  >  3  x  1012 

s'1.  These  enhanced  barrier-crossing  frequencies  are  roughly  consistent  with 
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some  theoretical  predictions.  They  are  also  in  harmony  with  some  very 
recent  information  on  real-time  polar  solvation  dynamics  extracted  from 
subpicosecond  fluorescence  measurements;  a  detailed  comparison  between  r#fI 
obtained  from  electron- transfer  and  time -resolved  fluorescence  measurements 
will  be  featured  in  a  forthcoming  article. 

A  related  illustration  of  the  variations  in  the  degree  of  electron¬ 
tunneling  versus  solvent -dynamical  control  as  the  solvent  is  altered  is 
provided  in  Table  III,  which  lists  ratios  of  k#x  values  for  each  metallocene 
versus  that  for  Cp2Co+/0,  k#x/k#x(Cp2Co) ,  in  seven  solvents.  The  advantage  of 
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this  format  is  that  no  knowledge  of  the  solvent  dependence  of  the  free-energy 
barriers  is  required,  since  this  should  largely  cancel  when  the  rate  ratios 
are  taken.  The  variations  in  these  ratios  with  the  solvent  therefore  provide 
a  monitor  of  the  extent  to  which  the  relative  barrier-crossing  frequencies  for 
the  different  reactions  are  sensitive  to  the  solvent  dynamics.  To  facilitate 
comparison,  the  solvents  are  again  listed  in  order  of  decreasing  r^1,  although 
data  in  PC  are  inserted  in  a  position  commensurate  with  the  value 
inferred  above,  2  x  1012  s'1.  The  increasing  emergence  of  electron- tunneling 
control  of  k#x  as  the  r^1  value  for  a  given  redox  couple  is  increased,  or  the 
redox  couple  is  altered  in  the  sequence  Cp2Co+/0  <  Cp*Co+/o  <  Cp2Co+/° 
<  CpjFe’*70  <  Cp2Fe+/°  <  HMFc+/0,  can  be  discerned  clearly  in  the  progressive 
decreases  in  the  k#x/k<x (Cp2Co)  ratios  observed  under  these  conditions  (i.e. 
up  the  columns  and  across  the  rows,  respectively,  in  Table  111). 

Extraction  of  Electronic  Coupling  Matrix  Elements 

A  relatively  direct  means  of  extracting  the  extent  of  electronic  coupling 
for  each  reaction  from  the  results  in  Fig.  1  involves  matching  them  with  a 
corresponding  sequence  of  curves  obtained  using  kjx  values,  k^x(calc), 
calculated  from  a  theoretical  model  for  a  suitable  sequence  of  H12  values. 
Figure  2  displays  such  a  set  of  log  k'x(calc)  -  log  curves  obtained  for 
the  same  solvents  as  in  Fig.  1.  A  detailed  presentation  of  the  underlying 
theory  forms  the  subject  of  ref.  lib.  (See  also  the  Appendix  below.) 
Briefly,  the  k,x(calc)  values  in  each  solvent  were  calculated  by  utilizing 
Eq.  (5)  with  the  spatial  integration  performed  from  a  distance  of  closest 
approach,  re,  of  7.6  A  (i.e.  twice  the  metallocene  radius,  3.8  A^a) .  The 
dependence  of  the  "local"  (unimolecular)  rate  constant,  k#t,  upon  the  reactant 
internuclear  separation,  r,  required  in  Eq.  (5),  was  obtained  as  follows.  The 
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primary  input  parameters  are  rL,  H12,  and  AG  .  While  rL  can  be  considered  to 

39  * 

be  r- independent  for  the  present  purposes,  both  H12  and  AG  are  sensitive  to 
the  reactant  pair  geometry. 

The  r  dependence  of  the  former  can  be  taken  as^ 

(H12)2  -  (H®2)2  expt-a  (r  -  r0)]  (7) 

where  Hj2  is  the  matrix  coupling  element  for  the  "closest  approach"  geometry. 

2  llh 

The  coefficient  a  is  taken  here  to  be  the  typical  value  1.25  A'1.  ’  The 

corresponding  r-dependent  barrier  height  is  obtained  from^^ 

AG*  -  AG*  -  H12  (8) 

In  addition  to  the  increases  in  AG  with  increasing  r  due  to  corresponding 
decreases  in  H12,  AGe  and  hence  AG  should  also  increase  under  these 
conditions  according  to^^ 

AG*(r)  -  AG*(ro)  [2a(a' 1  -  r*1)]  (9) 

where  a  is  the  reactant  radius  and  AG*(r0)  is  the  cusp  barrier  height  when 
r  -  r0  (-  2a) .  This  latter  quantity  was  set  equal  to  the  AG*  values  given  for 
each  solvent  in  Table  I.  Although  the  radial  dependence  of  AG*  in  Eq.  (9)  has 
its  origin  in  the  two-sphere  dielectric  continuum  treatment, ^  it  is  also 
approximately  consistent  with  r-dependent  E0J)  measurements  for  biferrocene 
cations^  and  related  systems. Strictly  speaking,  Eq.  (9)  presumes  that 
AGe  arises  entirely  from  solvent  reorganization,  i.e.  the  inner-shell 
(reactant  vibrational  distortional)  barrier,  AG*t,  is  small.  Although  our 
earlier  estimates  for  AG*g  for  these  metallocene  couples^ ’ ^  are  marred  by  a 
systematic  error  in  the  bond-distance  estimates,  it  is  likely  that 


AG*#  »  0.5  -  1  kcal  mol*1.47  (Note  that  the  Eop  estimates  of  AG*  necessarily 
contain  any  inner-shell  contribution  to  the  blnuclear  redox  system.) 

In  essence,  the  calculation  of  k#fc(r)  and  hence  k^z(calc)  from  these 
parameters  involves  evaluating  the  corresponding  adiabatic  rate  constant 
k*t(r),  and  the  corresponding  modifications  to  account  for  nonadiabatic 
barrier  crossing.  The  latter  was  accomplished  here  by  using  Eq.  (Al)  in  the 
Appendix;  note  that  generally  k,t  £  kjt.  These  nonadiabatic  effects  (k#1  <  1) 
are  felt  increasingly  for  smaller  H12  and/or  larger  r*1.  The  calculation  of 
k*t  takes  into  account  the  influences  of  barrier-top  roundedness  (i.e.  non¬ 
zero  H12)  as  well  as  the  effect  of  solvent  friction  (via  rL)  upon  the  barrier - 
crossing  frequency  i/n.  Equation  (Al)  also  includes  the  effects  of  solvent 
inertia;  these  yield  deviations  from  pure  "overdamped"  behavior  (where 
a  r'1)  anticipated  in  low-friction  media  where  i/n  approaches  the 
transition-state  theory  (TST)  limit,  whereupon  equals  the  solvent  inertial 

frequency  wo/2ir  and  k*t  -  kTST . 14  ’  25  (See  Appendix  and  ref.  lib.) 

The  six  sets  of  points  shown  in  Fig.  2,  with  "best  fit"  curves  drawn 
through  each  set,  correspond  to  the  H°2  values,  from  0.05  to  1.0  kcal  mol*1, 
as  indicated.  (The  particular  solvents  are  identified  by  the  same  numbering 
scheme  as  in  Table  I  and  Fig.  1).  By  and  large,  the  calculated  curves 
(Fig.  2)  show  a  striking  similarity  to  the  corresponding  experimental  plots  in 
Fig.  1,  in  that  the  log  k,x(calc)  -  log  r^1  slopes  decrease  for  smaller  H°2 
and/or  larger  r*1 ,  reflecting  increasing  reaction  nonadiabaticity .  Indeed,  it 
is  interesting  to  note  that  even  the  scatter  in  the  calculated  points  in 
Fig.  2  for  solvents  having  the  largest  r*1  values  (D20,  acetone,  acetonitrile, 
and  nitrome thane)  is  mimiced  reasonably  well  by  the  experimental  data  in 
Fig.  1.  (This  scatter  arises  from  the  onset  of  an  inertial  limit  upon  k#x 
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which  is  solvent  dependent,  wo  being  largest  for  D20  and  smallest  for 

„  25. 

acetone  ) . 

We  therefore  have  reasonable  confidence  in  extracting  approximate 
estimates  of  H°2  for  the  various  experimental  systems  by  matching  the  shapes 
of  the  experimental  and  calculated  log  k^x  -  log  i*1  curves,  together  with  the 
vertical  displacement  between  them.  This  procedure  yields  the  following  H°2 
estimates  (kcal  mol'1):  Cp2Co'f/o,  1.0;  Cp2Co4/0,  0. 5-1.0;  Cp2Co+/o,  0.25; 

Cp2Fe+/o,  0.2;  Cp2Fe+/0,  0.1;  HMFc+/°,  0.075.  As  a  pictorial  illustration  of 
these  differences,  Fig.  3  shows  a  pair  of  free  energy-reaction  coordinate 
profiles,  calculated  as  in  ref.  11a,  appropriate  for  AG*  -  5  kcal  mol'1  with 
H12  -  1.0  (solid  curve)  and  0.1  kcal  mol'1  (dashed  curve). 

The  above  H“2  estimates  for  the  cobaltocene  systems,  especially  Cp2Co+/0, 
may  well  be  significantly  smaller  than  the  actual  values  given  the  possibility 

that  AG*t  for  these  couples,  are  slightly  (ca  0.5  kcal  mol'1)  larger  than  for 
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the  ferrocene  systems.  The  presence  of  such  increased  AGlj  not  only  will 

decrease  k^  in  a  given  solvent,  but  also  may  diminish  significantly  the 

influence  of  solvent  friction,  i.e.  depress  the  log  k^x  -  log  r^1 
15  52 

dependence.  ’  At  least  the  latter  complication,  however,  should  disappear 
in  the  dynamically  most  rapid  solvents,  where  r^1  (or  u>0/2n)  approaches  the 
vibrational  frequencies,  ca  6  to  9  x  1012  s*1,^  characterizing  the  symmetric 
Cp-M-Cp  inner-shell  distortions.^""*  If  AG*t  is  indeed  0.5  kcal  mol*1  larger 
for  Cp2Co"f/0  than  for  Cp2Fe+/o,  then  to  correct  for  this  difference  the 
log  k^x  values  for  the  former  couple  will  need  to  be  increased  by  0.35  (i.e. 
ca  2.5  fold  in  kxx)  at  the  largest  log  r^1  values  in  Fig.  1.  This  would  yield 
a  log  k^x  difference  of  1.0  rather  than  ca  0.65  for  these  couples  in,  say, 
acetonitrile  and  acetone.  Assuming  that  H°2  *  0.1  kcal  mol"1  for  Cp2Fe+/o,  it 
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can  be  deduced  from  Fig.  2  that  this  correction  yields  an  approximately 
twofold  enhancement  in  the  H°2  estimate  for  Cp2Co+/°,  from  ca  0.25  to  0.5  kcal 
mol"1 . 

The  ca  3  fold  larger  kxx  values  for  Cp2Co+/0  relative  to  Cp2Co+/0 

observed  in  each  solvent  may  also  be  attributable  in  part  to  differences  in 
*  47 

AGlg.  This  behavior  is  also  consistent  with  a  smaller  outer-shell  (solvent 

reorganization)  barrier  for  the  former  couple  that  might  be  anticipated  from 

its  slightly  larger  effective  radius.  Alternatively,  H°2  for  Cp2Co+/o  may  be 

even  larger  than  1.0  kcal  rno'l*1.  Apart  from  the  uncertainties  in  extracting 

the  solvent -dependent  k^x  values  themselves,  the  estimated  values  are 
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subject  to  significant  error.  The  assumption  that  the  first  eleven 

solvents  in  Table  1  exhibit  Debye  behavior  is,  of  course,  only  an 

approximation,  even  though  the  most  recent  subpicosecond  time -resolved 

fluorescence  data  yields  major  relaxation  times  that  typically  differ  from  rL 
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by  twofold  or  less.  [The  observed  positive  deviations  of  the  k^x  values 

from  the  anticipated  log  k^x  -  log  r^1  plots  in  benzonitrile ,  for  example 
(Fig.  1),  may  be  due  to  the  influence  of  a  higher -frequency  relaxation  in  this 
solvent.]  Possible  solvent-specific  effects  upon  the  approach  of  the  reacting 
partners,  which  would  be  reflected  in  non-unit  gDA(r)  values  [Eq.  (5)],  may 
also  contribute  significantly  to  the  rate-solvent  dependence.  Generally 
speaking,  however,  consideration  of  these  and  related  factors  lead  only  to 
relatively  minor  (<  2  fold)  uncertainties  in  the  inferred  H°2  estimates.  The 
relative  HJ2  values,  then,  are  deemed  sufficiently  reliable  to  warrant  at 
least  semiquantitative  interpretation. 

Interpretation  and  Implications  of  H°:  Values 

The  markedly  (ca  2.5-5  fold)  larger  H°2  values  estimated  for  Cp2Co+/° 
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and  other  cobaltocene  couples  compared  with  their  ferrocene  analogs  can  be 
accounted  for  in  terms  of  the  spatial  differences  in  the  "redox  molecular 

C 

orbitals”,  corresponding  to  the  UIMO/HOMO  pair  for  the  Cp2M+/cP2M  reaction 

partners.^  While  theoretical  studies  indicate  that  the  cobaltocene  HOMO  is 
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delocalized  over  the  Cp  rings,  that  for  ferrocene  is  strongly  metal 

centered. ^  Since  similar  considerations  also  apply  to  the  LUMO  for  the 
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oxidized  redox  forms,  '  *  the  donor-acceptor  orbital  overlap  should  be 


greater 

for 

the 

Cp2Co+/0  versus 

the 

Cp2Fe+/0  couples.^  The 

relative  H°z 

values 

for 

these 

reactions  are 

in 

harmony  with  the  twofold 

greater  H12 

estimates  derived  from  optical  electron  transfer  data  for  bicobaltocene  versus 

ft 

biferrocene  cations,  although  the  absolute  values  of  the  latter  are  much 
larger  (3.2  and  1.7  kcal  mol-1,  respectively)  partly  as  a  result  of  direct 
bond  formation  between  the  redox  centers. 

Very  recent  ab  initio  calculations  by  Newton  also  yield  markedly  (ca  3-6 
fold)  larger  HJ2  values  for  Cp2Co+/0  than  for  Cp2Fe+/0  self  exchange.®  These 
calculations  also  indicate  that  considerably  stronger  overlap  is  obtained  for 
approach  of  the  reacting  pair  along  a  common  fivefold  axis  ("axial"  geometry) 
than  for  a  "side-by-side"  configuration.  For  the  former  geometry,  H12  -  2.5 
and  0.4  kcal  mol-1  for  Cp2Co+/0  and  Cp2Fe+/0,  respectively,  when  the  Cp  rings 
are  "in  contact"  (3.5  A  apart);  whereas  in  the  latter  configuration, 
H“2  «  0.04  and  0.1  kcal  mol-1,  respectively.  These  findings  foreshadow  a 
limitation  in  the  present  analysis;  the  substantial  anisotropy  in  the 
calculated  H12  values  suggests  'that  the  assumption  of  spherical  symmetry 
embodied  in  Eq.  (5)  is,  strictly  speaking,  inadequate.  Nevertheless,  the 
present  experimental  H°2  estimates  should  represent  at  least  weighted  averages 
of  the  true  values  for  the  various  "close  contact"  encounter  pair  geometries, 
and  as  such  are  reasonably  consistent  with  the  calculated  values.  It  would  be 
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of  interest  to  perform  more  sophisticated  kinetic  analyses  employing  such 
detailed  theoretical  information. 

Interpretation  of  the  substantial  observed  variation  in  H°2  resulting 

from  Cp  ring  substitution,  unlike  alteration  of  the  metal  center,  is  limited 

by  the  lack  of  electronic  structural  information  for  these  systems.  At  first 

sight,  one  night  expect  that  H°2  would  correlate  roughly  with  empirical 

electronic  quantities  such  as  the  Hammett  a  parameter  which  characterizes 
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substituent  effects  in  organic  chemistry,  especially  since  they  correlate 

well  with  metallocene  formal  potentials^’ ^  and  NMR  chemical  shifts. ^ 

To  this  end,  Table  IV  contains  ktx  values  in  acetonitrile  for  the  self- 

exchange  reactions  studied  here,  together  with  some  additional  data  for  four 

monosubstituted  ferrocenes , ^  in  comparison  with  corresponding  formal 

potentials,  E£.  (The  entries  are  listed  in  order  of  increasing  Ef.)  As  noted 

above,  the  differences  in  k#x  in  this  dynamically  rapid  solvent  should  be 

particularly  sensitive  to  variations  in  Hl2.  The  lack  of  any  correlation 

between  k#x  and  E£  for  either  the  cobaltocene  and  ferrocene  couples  is  clearly 

evident,  even  though  substantial  variations  in  both  these  quantities  are 

observed.  Upon  reflection,  however,  this  finding  is  unsurprising  since  the 

substituent  effect  upon  E£  arises  from  alterations  in  the  energy  of  the  redox 

orbital,  whereas  Hl2  and  hence  k>x  is  responsive  primarily  to  the  spatial 
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properties  of  the  HOMO/LUMO  pair.  Given  the  magnitude  of  such  apparent 
substituent  effects  upon  H12,  their  further  detailed  examination  would  be  of 
great  interest. 

As  noted  previously,^  the  sensitivity  of  the  self -exchange  kinetics  to 
the  metallocene  structure  contrasts  with  the  exchange  kinetics  of  these 
systems  at  metal  electrodes.  In  the  latter  environments,  uniformly  facile 
kinetics  are  obtained  that  apparently  exhibit  a  strong  dependence  upon  the 
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solvent  dynamics,  Indicative  of  essentially  adiabatic  behavior . ^ ^  The 

likelihood  that  this  behavior  persists  for  ferrocenium- ferrocene  couples  at 

electrodes,  even  in  solvents  such  as  acetonitrile, ^  suggests  that 

electrochemical  exchange  processes  may  proceed  by  inherently  more  adiabatic 

pathways  than  for  related  self -exchange  reactions.  Indeed,  such  a  conclusion 
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is  in  harmony  with  the  prediction  of  recent  theoretical  treatments . 

Although  detailed  theoretical  calculations  of  H°z  have  only  been 
undertaken  so  far  for  a  relatively  narrow  range  of  homogeneous -phase  self- 
exchange  reactions,  most  values  are  substantially  below  both  the  magnitude, 
«c  0.2  kcal  mol'1,  necessary  to  maintain  a  significant  degree  of  reaction 
adiabaticity  in  dynamically  rapid  solvents,  such  as  water.  At  least  in 
homogeneous  solution,  then,  the  range  of  outer- sphere  reactions  influenced 
importantly  by  solvent  friction  effects  may  well  be  smaller  than  those  subject 
primarily  to  electron- tunneling  control  of  the  barrier-crossing  frequency. 
Although  the  present  solvent -dependent  analysis  is  applicable  only  in 
unusually  favorable  circumstances,  it  would  be  of  considerable  interest  to 
utilize  it  to  explore  electronic  coupling  factors  in  other  types  of  processes, 
especially  involving  intramolecular  electron  transfer. 
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Appendix 

A  key  element  in  obtaining  the  calculated  k^x  values,  k^x(calc),  shown  in 
Fig.  2  involves  the  computation  of  unimolecular  rate  constants,  k>t(r),  that 
span  the  full  range  from  entirely  nonadiabatic  to  adiabatic  behavior.  For 
this  purpose,  we  use  the  following  relation,  derived  in  ref.  lib: 
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where  kjt(r)  is  the  rate  constant  for  adiabatic  barrier  crossing,  and  kTST  is 
the  rate  constant  corresponding  to  the  TST  limit  (i.e.  in  the  absence  of 
solvent  friction)  given  by 

kTST  -  <w0/2ir)  exp(-AG*/RT)  A2) 


where  «o  is  the  solvent  inertial  frequency.  The  electronic  transmission 
coefficient,  <t#1,  can  be  related  to  the  more  familiar  <c#1  in  Eq.  (1)  by  (see 
Appendix  to  ref.  lib) 


*.i  “  +  “.id  -  ".)]  (A3) 

where  the  "adiabatic  transmission  coefficient"  ica  -  k*t/kTST- 

Apart  from  the  additional  -Jc#1  term  in  the  denominator,  Eq.  (Al)  is 
similar  or  identical  to  other  literature  relations  for  interpolating  between 
the  adiabatic  and  nonadiabatic  limits  (e.g.  refs.  12,  13).  This  term  allows 
k#t(r)  to  reduce  to  the  correct  adiabatic  TST  limit,  whereupon  the 
preexponential  factor  equals  w0/2ir.^^ 

The  adiabatic  rate  constant  kjt  was  calculated  using  a  position  (q)- 
dependent  diffusional  model,  whereby  the  diffusion  coefficient,  D,  describing 
progress  along  the  reaction  coordinate  in  the  presence  of  solvent  friction 
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within  the  intersection  region  differs  from  that  in  the  potential -energy  wells 
as  a  result  of  barrier-top  roundedness  (i.e.  H12  >  0).*^  Specifically,  k*t 
was  approximated  by 

<Kt r1  -  2  /  dq  (D(q)  P.^q)]"1  (A4) 

O 

where  the  integration  is  performed  from  the  reactant  well  to  the  barrier  top; 
P.q(q)  is  given  by 

P.qCq)  *  (2*//J  w2)*1  exp[  -f)  V(q>]  (A5) 

where  fl  -  (kjT)'1  (kj  is  the  Boltzmann  constant)  and  V(q)  describes  the 

potential -energy  surface.  Appropriate  expressions  for  D(q)  in  terms  of  rL 

(i.e.  for  Debye  solvents)  and  for  V(q)  are  given  in  ref.  lib,  along  with  other 

pertinent  details.  The  effect  of  increasing  H12  is  to  decrease  k*t;  typically 

k*t  decreases  by  2-3  fold  as  Hl2  increases  from  0  to  1  kcal  mol"1,  the  effect 

depending  somewhat  on  the  extent  of  solvent  friction.  The  results  are 
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numerically  similar  to  those  obtained  by  utilizing  the  Smoluchowski 

14 

treatment  of  Calef  and  Wolynes.  (Note  that  the  effect  of  solvent  friction 
upon  k-t,  as  given  by  Eq.  (Al),  is  contained  entirely  within  the  kjt  term.) 

The  JE#1  values  were  calculated  by  using  the  approximate  expression'^ 


lc9l  -  4w3/2  7r/(l  +  7r) 

(A6) 

where 

7t  -  |H12|2/4A  *0(AG*  kjT)** 

(A7) 

For  the  present  purposes,  it  is  most  important  to  note  that  the  derived  values 
of  and  hence  the  form  of  the  log  k>x(calc)  -  log  r"1  curves  in  Fig.  2, 
are  insensitive  to  the  particular  expression  used  for 
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T.  Summary  of  Rate  Constants,  k#Jt,  for  Metallocene  Self  Exchange  Reactions  at  25*C,  and 
Related  Solvent-Dependent  Data 

10'7  kM,  M'1  s'1  d 


Solvent* 

1012  s*1 

AG^  Cp2Co+/°e  CP;Co+/°f 

Cp2Co*/og 

Cp^Fe+/oh  Cp2Fe+/oh 

HMFc+/of 

1.  Acetonitrile 

5.35 

55 

20 

4.5 

3.5 

0.9 

0.55 

2.  Propionitrile 

_4-M 

5.15 

5.5 

3.1 

0.92 

-0.9 

3.  Acetone 

3. 5k 

5.4 

28 

12 

2.3 

2.2 

0.8 

0.72 

4.  D20 

1.91 

5.2 

«  30 

1.4 

5.  Nitrome thane 

4.5“ 

5.3 

24 

3.7 

1.2 

-1.4 

6.  DMF 

0 . 77n 

(5.0) 

25 

3.2 

7 .  DMSO 

0.5° 

(4.9) 

23 

3.0 

0.95 

8.  Benzonitrile 

0.2P 

4.55 

30 

8 

7.5 

6 

2.7 

1.8 

9.  Nitrobenzene 

0.2P 

4.55 

S  4.5 

6 

3.0 

1.4 

10.  TMU 

0.179 

(4.85) 

4.9 

1.55 

11.  HMPA 

0.11n 

(4.95) 

4.3 

0.67 

12.  PC 

(0.4)r 

5.25 

35 

14 

3.6 

1.2 

13.  NMF 

0 . 27n 

(5.15) 

8.5 

14.  Methanol 

(0.135) 

s  5.2 

-15 

9 

1.8 

1.7 

15.  Ethanol 

(0.033) 

t  5.0 

6 
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a  DMF  -  N,N- dime thylformami.de ,  DMSO  -  dimethylsulfoxide ,  PC  -  propylene 
carbonate,  NMF  -  N-methylformamide ,  TMU  -  tetramethylurea,  HMPA  - 
hexamethylphosphoramide . 

^  Inverse  of  longitudinal  solvent  relaxation  time,  obtained  from  Eq.  (3) 
using  values  of  rD  and  fa  given  in  literature  values  cited,  and  c0  either 
from  these  sources  or  ref.  23.  (In  most  cases,  these  values  coincide  with 
tabulations  in  refs.  5a  and  6c.  Vhere  dielectric  values  not  listed  in 
refs.  5a  or  6c  are  employed  to  estimate  rL,  these  are  given  following 
appropriate  footnote  to  literature  citation.) 

c  Free-energy  barrier  (in  "cusp  limit",  H12  -*  0)  for  self-exchange 
reactions  in  a  given  solvent,  estimated  from  energy  of  intervalence  band 
maximum,  E  ,  in  near-infrared  spectra  for  biferrocenylacetylene  cation 
using  Eq.  4.  Value  for  D20  estimated  from  corresponding  E  value  for 
biferrocene  cation  (see  ref.  5c).  Values  in  parentheses  were  estimated  as 
described  in  the  text. 

^  Self- exchange  rate  constant  of  various  cobalt  and  ion  metallocene  redox 
couples  in  solvent  indicated,  as  obtained  by  proton  NMR  line  broadening. 
Cp  -  cyclopentadienyl ,  Cp'  -  pentamethylcyclopentadienyl ,  Cp*  - 
monocarboxymethylcyclopentadienyl ,  HMFc  -  hexamethylferrocene .  All  values 
measured  at  (or  extrapolated  to)  total  ionic  strengths  n  «  0.01-0.02  {1. 
(See  ref.  5b  for  details  of  ionic -strength  dependence  of  k  .) 
Reproducibility  of  k#x  values  typically  within  ca  5-10%,  except  *for 
CpJCo’f/0  for  which  reproducibility  was  mostly  within  ca  20%  (see  text)  . 

e  Values  taken  from  data  in  ref.  5a,  but  increased  by  ca  25%  so  to  account 
for  decrease  of  n  from  0.14  to  0.02  JJ. 

^  Obtained  in  present  study  (also  ref.  5c). 

&  As  for  footnote  e,  but  k#x  values  in  propionitrile  and  ethanol  obtained 
in  present  study. 

^  Obtained  in  present  study,  and  from  ref.  5b. 

*  Estimated  from  data  in  ref.  24. 

J  Effective  rL  value,  estimated  by  accounting  for  minor  non-Debve 
characteristics  of  dielectric  loss  spectra,  appearing  as  additional  loss 
components  at  frequencies  above  rj1  (also  see  ref.  25). 

If 

Estimated  from  data  in  ref.  26  (also  see  ref.  25). 

1  Ref.  27  ®  Ref.  28  n  Ref.  29 

°  Ref-  30  (td  -  19  ps,  -  5.1,  e0  -  46.7). 

p  Ref.  31  9  Ref.  32  r  Ref.  33  s  Refs.  34,  35  c  Ref.  34 
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TABLE  II.  "Barrier-Corrected"  Rate  Constants  for  Metallocene  Self  Exchange  in 
Various  Solvents  Relative  to  those  in  Acetonitrile,  k^x/k^x(ACN) ,  in 
Comparison  with  Corresponding  Relative  Solvent  Relaxation  Dynamics, 
r’Vr-^ACN) 


Solvent- 

r^/f^CACN)6 

k;xA;x 

(ACN)C 

Cp2Co+/0 

Cp;Co+/0 

Cp2Co+/° 

Cp2Fe+/0 

Cp2Fe+/° 

HMFc+/° 

Propionitrile 

-  1 

0.9 

0.65 

0.75 

1.2 

Acetone 

0.9 

0.6 

0.7 

0.6 

0.7 

1.0 

1.4 

d2o 

0.5 

1.2 

2 

Nitrome thane 

*  1.1 

1.1 

1.0 

1.25 

2.3 

DMF 

0.2 

0.25 

0.13 

0.4 

DMSO 

0.12 

0.2 

0.3 

0.5 

Benzonitrile 

0.05 

0.15 

0.11 

0.4 

0.45 

0.8 

0.9 

Nitrobenzene 

0.05 

0.06 

0.45 

0.9 

0.7 

TMU 

0.043 

0.11 

0.15 

HMPA 

0.028 

0.04 

0.075 

PC 

(0.1) 

0.55 

0.6 

0.65 

1.9 

NMF 

0.07 

1.4 

Methanol 

(0.035) 

0.6 

1.6 

1.55 

2.4 

Ethanol 

(0.008) 

0.75 

a  See  footnote  a  to  Table  I  for  explanation  of  abbreviations. 

®  Ratio  of  inverse  longitudinal  relaxation  time  for  given  solvent  to  that 
in  acetonitrile.  See  Table  I  for  r^1  values  and  data  sources. 

c  Ratio  of  barrier-corrected  rate  constant,  k^x,  for  given  metallocene 
redox  couple  in  various  solvents  to  that  in  acetonitrile,  k^x(ACN).  Values 
of  k^x  obtained  from  experimental  k#x  and  AG*  values  (Table  I)  by  using 
Eq.  (6).  See  footnote  d  of  Table  I  for  explanation  of  metallocene 
abbreviations . 
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TABLE  III .  Ratios  of  Rate  Constants  for  Metallocene  Self  Exchange  to 
those  for  Cp2Co+/°,  k<x/ktJt(Cp2Co)  ,  in  Various  Solvents 


Solvent3 

10‘V‘ 

A„(CpiCo)c_ 

Cp2Co+/° 

Cp2Co+/° 

Cp2Fe+/0 

Cp2Fe+/o 

HMFc+'° 

Acetonitrile 

4 

0.35 

0.08 

0.065 

0.015 

0.010 

Acetone 

3.5 

0.43 

0.08 

0.08 

0.03 

0.025 

PC 

(0.4) 

0.4 

0.11 

0.035 

DMF 

0.77 

0.13 

DMSO 

0.5 

0.13 

Benzonitrile 

0.2 

0.27 

0.25 

0.20 

0.09 

0.06 

HMPA 

0.11 

0.16 

a  See  footnote  a  to  Table  I  for  explanation  of  abbreviations, 
k  Inverse  longitudinal  solvent  relaxation  time  (as  in  Table  I). 
c  Calculated  from  k#3t  values  in  Table  I . 
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tabt.f  tv.  Comparison  between  Self -Exchange 
Metallocene  Redox  Couples  in 
Corresponding  Formal  Potentials 


Rate  Constants  for  Various 
Acetonitrile  at  25°C  and 


Redox  Couple* 

Et b 

V.  vs.  Cp2Fe+/0 

k  C 

•X 

H'1  s*1 

Cp2Co+/° 

-1.92 

5.5  x  10® 

Cp2Co+/° 

-1.30 

4.5  x  107 

Cp*Co+/0 

-0.79 

2.0  x  10® 

Cp2Fe+/° 

-0.50 

3.5  x  107 

(Cp-CpCH20H)Fe+/o 

-0.01 

5.5  x  106 

Cp2Fe+/0 

0 

9  x  106 

(Cp  ‘  CpI)Fe+/° 

0.15 

9.2  x  106 

(Cp-CpC0CH3)Fe+/o 

0.24 

8.0  x  106 

(Cp • CpCN)Fe+/° 

0.36 

5.7  x  10® 

a  Nomenclature  partly  as  in  Tables  I -III.  [Note  that 
hydroxymethyl (ferrocenium- ferrocene)  is  abbreviated  as  (Cp-Cp‘CH20H)Fe+/o 
rather  than  HMFc+/°  for  consistency  with  other  monosubstituted  ferrocenes 
listed  here . ] 

®  Formal  potential  (volts)  versus  ferrocenium- ferrocene  in  same  solvent 
(acetonitrile)  evaluated  at  gold  or  mercury  electrodes  in  0.1  & 
tetrabutylammonium  hexafluorophosphate  by  using  cyclic  voltammetry. 

c  Rate  constant  for  metallocene  self -exchange  in  acetonitrile,  for  ionic 
strength  n  -  0-0.02  Hi  taken  from  Table  I  and  ref.  57.  Values  reproducible 
generally  to  ±  10-20%. 
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Logarithmic  plots  of  "barrier -corrected"  rate  constants,  kxx,  versus 
inverse  of  longitudinal  relaxation  time ,  ,  for  five  metallocene  redox 
couples  in  eleven  solvents.  kxx  values  were  extracted  from  experimental 
k#x  and  AG*  values  given  in  Table  1  by  using  Eq.  (6).  See  Table  1  for 
solvent  numbering  scheme.  Key  to  redox  couples:  Cp^Co-*70,  filled  circles; 
CpjCo+/°,  filled  squares;  Cp  2Co+/0,  filled  triangles;  Cp2Fe+/°,  open 
triangles;  HMFc+/0,  open  squares.  Solid  traces  through  each  set  of 
solvent -dependent  log  kxx  values  are  suggested  "best  fit"  curves  (see 
text) . 

Fig.  2 

Logarithmic  plots  of  calculated  "barrier-corrected"  rate  constants, 
kxx(calc),  versus  inverse  of  longitudinal  relaxation  time,  r'1,  in  same 
eleven  solvents  as  for  Fig.  1,  for  the  sequence  of  six  electronic  matrix 
coupling  elements  at  reactant  "contact",  H°2,  as  indicated.  Values  of 
k^x(calc)  obtained  from  corresponding  calculated  values,  ktx(calc),  in  each 
solvent  by  using  Eq.  (6)  in  the  same  manner  as  the  experimental  quantities. 
Values  of  k#x(calc)  obtained  by  using  spatial  integration  procedure  [Eq. 
(5)]  with  rQ  taken  as  7.6  A,  and  gDA(r)  -  1.  Required  values  of  k#t(r) 
obtained  from  Eq.  (Al)  with  constituent  quantities  calculated  as  follows: 
kIST  from  Eq.  (A2)  using  "solvent  inertial"  frequencies,  u0,  calculated  as 
described  in  ref.  25  (see  footnotes  to  Fig.  4  and  Table  II  of  ref.  25  for 
numerical  values);  kjt(r)  from  Eqs.  (9) - (16)  of  ref.  lib;  «#1  from  Eqs. 
(20),  (20a)  of  ref.  lib.  The  r-dependent  AG*  values  required  for  kISI  and 
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k*t(r)  were  obtained  from  Eqs.  (8)  and  (9)  with  AG*(ro)  as  given  in  Table  I 

and  a  -  3.8  A.  The  r-dependent  H12  values  also  required  in  Eq.  (8),  as 

well  as  for  calculating  ic#l,  were  obtained  from  Eq.  (7)  with  a  -  1.25  A*1. 

The  double  integrals  required  using  this  overall  procedure  (both  spatially 

and  along  the  reaction  coordinated^)  were  evaluated  using  a  standard  two 

64 

dimensional  quadrature  technique,  executed  on  the  Purdue  University  Cyber 
205  supercomputer  using  double  precision  (128  bit). 

ru,-3 

Schematic  free  energy  -  reaction  coordinate  profiles  for  electron- 
exchange  reactions  having  different  electronic  matrix  coupling  elements, 
H12.  Curves  calculated  by  using  Eqs.  (4)  and  (5)  of  ref.  11a.  For  a 
5  kcal  mol'1  cusp  barrier,  the  solid  and  dashed  traces  shown  correspond  to 
H12  values  of  1.0  and  0.1  kcal  mol"1,  respectively. 
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NMR  Parameters  and  Related  Data  for  Metallocene  Self  Exchange 


Solvent3 

[Cp2M+]c 

fCp2M]d 

W  e 

WDP 

Wf 

wp 

c  h 

6d 

Sp1 

O 

7? 

• 

H 

mil 

mil 

Hz 

Hz 

ppm 

ppm 

ppm 

tr1  s'1 

HMFc+/o 

Self  Exchange  at 

200.0  MHz 

and  25° 

c  k 

Acetonitrile*3 

1.4 

36.6 

63 

287 

3.41 

4.30 

-20.18 

0.54 

Propionitrile 

13.5 

34.6 

232 

319 

-2.78 

4.27 

-20.85 

0.94 

Acetone*3 

1.1 

33.8 

45.7 

317 

3.59 

4.33 

-20.27 

0.72 

d2o 

4.5 

5.4 

751 

394 

-7.29 

4.33 

-21.33 

1.45 

Nitrome thane 

11.35 

31.2 

179 

260 

-2.34 

4.47 

-21.04 

1.4 

Benzonitrile 

5.1 

33.7 

102 

403 

1.39 

4.54 

-20.05 

1.8 

Nitrobenzene*3 

0.95 

26.1 

44.3 

450 

3.68 

4.54 

-20.47 

1.4 

Methanol*3 

7.5 

27.4 

177 

360 

-1.30 

4.35 

-21.21 

1.7 

Cp2Co+/0 

Self  Exchange  at 

469.5  MHz 

and  25 

“C  b 

Acetonitrile*3 

23.9 

14.2 

27 

38 

8.30 

3.91 

15.81 

19 

Acetone*3 

22.6 

10.6 

31 

32 

7.49 

3.83 

15.23 

12 

d2o 

23.1 

0.5 

3.5 

45 

4.17 

3.88 

16.80 

=30 

Nitrome thane 

17.3 

12.4 

32 

44 

8.89 

3.93 

15.80 

24 

Benzonitrile 

18.3 

2.1 

24 

40 

5.15 

3.95 

15.71 

9 

Nitrobenzene*3 

32 

0.75 

2.6 

51 

4.28 

4.01 

15.76 

-4.5 

TMU 

19.6 

10.0 

79 

63 

7.83 

3.95 

15.47 

4.9 

PC 

17.2 

15.8 

33 

60 

9.64 

3.93 

15.84 

14 

Methanol*3 

27.1 

11.5 

37 

54 

7.78 

3.96 

16.82 

12 

[See  next  page  for  footnotes] 


Footnotes 


a  TMU  -  tetramethylurea ,  PC  -  propylene  carbonate. 

^  Deuterated  solvent  used. 

c  Concentration  of  oxidized  form  of  redox  couple. 

^  Concentration  of  reduced  form  of  redox  couple . 

e  Linewidth  (at  half  height)  for  diamagnetic -paramagnetic  mixture. 

^  Linewidth  for  pure  paramagnetic  species  (linewidth  for  pure  diamagnetic 
species  taken,  or  measured  to  be,  1  Hz  in  most  cases). 

8  Chemical  shift  (ppm  vs.  TMS)  for  diamagnetic -paramagnetic  mixture. 

^  Chemical  shift  for  pure  diamagnetic  species. 

1  Chemical  shift  for  pure  paramagnetic  species. 

J  Rate  constant  for  self  exchange,  obtained  from  listed  NMR  parameters  as 
outlined  in  refs.  5a  and  b. 

HMFc  -  hydroxymethylferrocene ;  NMR  parameters  listed  refer  to  methylene 
proton  resonances  (see  Experimental  Section) . 

1  Cp'Co  -  carboxymethylcobaltocene ;  NMR  parameters  listed  refer  to  methyl 
proton  resonances  (see  Experimental  Section). 


